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A corpuscular simulation model of optical phenomena that does not require the knowledge of the solution of 
a wave equation of the whole system and reproduces the results of Maxwell's theory by generating detection 
events one-by-one is discussed. The event-based corpuscular model gives a unified description of multiple-beam 
fringes of a plane parallel plate and single-photon Mach-Zehnder interferometer, Wheeler's delayed choice, pho- 
ton tunneling, quantum eraser, two-beam interference, Einstein-Podolsky-Rosen-Bohm and Hanbury Brown- 
Twiss experiments. The approach is illustrated by application to a recent proposal for a quantum-controlled 
delayed choice experiment, demonstrating that also this thought experiment can be understood in terms of par- 
ticle processes only. 
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INTRODUCTION 

Quantum theory has proven extraordinarily powerful for 
describing the statistical properties of a vast number of labora- 
tory experiments. Conceptually, it is straightforward to use the 
quantum theoretical formalism to calculate numbers that can 
be compared with experimental data, at least if these numbers 
refer to statistical averages. However, a fundamental prob- 
lem appears if an experiment provides access to the individ- 
ual events that collectively build the statistical average. Prime 
examples are the single-electron two-slit experiment 01 > neu- 
tron interferometry experiments [2] and similar experiments 
in optics where the click of the detector is identified with the 
arrival of a single photon [3]. Although quantum theory pro- 
vides a recipe to compute the frequencies for observing events 
it does not account for the observation of the individual detec- 
tion events themselves JH|^]. For a recent review of various 
approaches to the quantum measurement problem and an ex- 
planation of it within the statistical interpretation, see Ref.0- 

From the viewpoint of quantum theory, the central issue 
is how it can be that experiments yield definite answers. As 
stated by Leggett [7]: "In the final analysis, physics cannot 
forever refuse to give an account of how it is that we obtain 
definite results whenever we do a particular measurement". 

This paper is not about interpretations or extensions of 
quantum theory. It gives a brief account of a very different 
approach to deal with the fact that experiments yield definite 
results. The latter, which is intimately linked to human per- 
ception, is taken as fundamental. We call these definite re- 
sults "events". Instead of trying to fit the existence of these 



events in some formal, mathematical theory, we change the 
paradigm by directly searching for the rules that transform 
events into other events and, by repeated application, yield fre- 
quency distributions of events that agree with those predicted 
by quantum theory. Obviously, such rules cannot be derived 
from quantum theory or, as a matter of fact, of any theory that 
is probabilistic in nature simply because these theories do not 
entail a procedure (= algorithm) to produce events themselves. 

The event-based approach has successfully been used to 
perform discrete-event simulations of the single beam split- 
ter and Mach-Zehnder interferometer experiment of Grang- 
ier et al. JUl (see Refs. jj- IT 1), Wheeler's delayed choice 



experiment of Jacques et al. [12] (see Refs. |11 




the quantum eraser experiment of Schwindt et aJ. 
Ref. Ill, [HI]), double-slit and two-beam single-photon in- 



terference experiments and the single-photon interference 
experiment with a Fresnel biprism of Jacques et aJ. lfl7ll 
(see Ref. Ill ll Il8l0 . quantum cryptography protocols (see 
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Ref. 1 19]), the Hanbury Brown-Twiss experiment of Aga- 
fonov et al. (see Ref. UHI^), universal quantum com- 
putation (see Ref. |22, 23|), Einstein-Podolsky-Rosen-Bohm- 
type of experiments of Aspect et al. B24 12511 and Weihs et 
al. 12611 (see Refs. Illll27l432ll ). and the propagation of elec- 
tromagnetic plane waves through homogeneous thin films and 
stratified media (see Ref. llllLl33ll ). An extensive review of the 
simulation method and its applications is given in Ref. [11]. 

A detailed discussion of the discrete-event approach can- 
not be fitted in this short paper. Therefore, we have chosen 
to illustrate the approach by an application to a recent pro- 
posal for a quantum-controlled Wheeler delayed choice ex- 
periment I134tl . We demonstrate that also this thought experi- 
ment can be understood in terms of event-based, particle-like 
processes only. The presentation is sufficiently detailed such 
that the reader who is interested can reproduce our results. 
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FIG. 1: Diagram of a standard Wheeler delayed-choice experiment 
with a Mach-Zehnder interferometer. Photons enter the interferom- 
eter via 50-50 beam splitter 1 (BS1). In the wave picture, the par- 
tial wave traveling along path (1) acquires a phase shift tj>o (<j>i). 
The variable x = 0, 1 controls the presence of 50-50 beam splitter 2 
(BS2). If BS2 is not in place (x = 0, indicated by the dashed rectan- 
gle) the partial waves do not interfere and the probability to observe 
the photon in path or 1 does not depend on the phase shifts. If 
BS2 is in place (x = 1, indicated by solid rectangle) the partial waves 
interfere and the probability to observe the photon in path or 1 is 
given by (1 +cos(0o — 0l))/2 or (1 — cos(0 o — 0i))/2, respectively. 
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FIG. 2: Quantum gate representation of the standard Wheeler 
delayed-choice experiment with a Mach-Zehnder interferometer (see 
Fig.[T}. The first Hadamard gate H acts as a 50-50 beam splitter, 
changing the state |0) into the state (|0) + 1 \))/sqrt2. The phase gate 
<p changes the amplitude of the state |1) by e'f. The second (con- 
trolled) Hadamard gate H act as a 50-50 beam splitter if the control 
variable x = 1 or passes the photons unaltered if x = 0. The angle a 
determines the probability that the control variable x is 1. A pair of 
detectors (not shown) signals the presence of a photon in the state |0) 
or 1 1) and with each detected photon the value of x is being recorded. 
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WHEELER'S DELAYED-CHOICE EXPERIMENT 



Particle-wave duality, a concept of quantum theory, at- 
tributes to photons the properties of both wave and particle 
behavior depending upon the circumstances of the experi- 
ment ||4j]. The particle behavior of photons has been shown 
in an experiment composed of a single beam splitter (BS) 
and a source emitting single photons and pairs of photons JUl- 
The wave character has been demonstrated in a single-photon 
Mach-Zehnder interferometer (MZI) experiment |8|]. The lay- 
out of such an experiment is shown in Fig. Q] By adding 
a device which controls the presence or absence of the sec- 
ond beam splitter BS2, this setup can be used to perform a 
delayed-choice experiment. Originally, Wheeler proposed a 
double-slit gedanken experiment in which the decision to ob- 
serve wave or particle behavior is made after the photon has 
passed the slits [35]. Similarly, in the MZI experiment, the 
decision to remove and place BS2 at the intersection of paths 
and 1 can, in principle, be made after the photon has passed 
BS1. The conclusion is that the pictorial description of this 
experiment defies common sense: The behavior of the photon 
in the past is said to be changing from a particle to a wave or 
vice versa. 



FIG. 3: Quantum gate representation of the quantum version of 
Wheeler delayed-choice experiment with a Mach-Zehnder interfer- 
ometer [34]. Reading from left to right, the first Hadamard gate H 
on the top line acts as a 50-50 beam splitter and the phase gate (p 
changes the amplitude of the state 1 1) by e'f. The second (controlled) 
Hadamard gate H on the top line acts as a 50-50 beam splitter if the 
state of the ancilla is 1 1) or passes the photons unaltered if that state 
is |0). Initially in the state |0), the ancilla is prepared in a uniform 
superposition of the states |0) and 1 1) by another interferometer cir- 
cuit (bottom line) in which the phase gate a changes the amplitude 
of the ancilla state 1 1) by e ia . The angle a determines the probabili- 
ties of the states |0) and 1 1). A pair of detectors (not shown) signals 
the presence of the photon in the state |0) or 1 1). Similarly, another 
pair of detectors (not shown) signals the presence of the ancilla in the 
state |0) or |1). 



QUANTUM CONTROLLED DELAYED-CHOICE 
EXPERIMENT 



It is of interest to enquire what happens if the variable x 
which controls the presence of BS2 (see Fig. [TJ or, equiva- 
lently, the controlled Hadamard gate (see Fig. |2]l is replaced 
by a quantum two-state system d 3411 . In a sense, one could 
then view the experiment as a simple example of a quantum- 
controlled experiment 113411 . The original proposal of the 
quantum-controlled delayed-choice experiment 13411 is formu- 
lated in a notation that is commonly used in the quantum com- 
puter literature ll3r3l . To facilitate the comparison with this 
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work, we also adopt to this notation from now on. First, in 
Fig. |2] we show the quantum gate diagram that is equivalent 
to the standard delayed-choice experiment depicted in Fig.Q] 
The main change, irrelevant from a conceptual point of view, 
is to replace the beam splitters by Hadamard gates. In Ref.b4 
it is proposed to replace the classical random variable x in 
Fig.|2]by a qubit, conventionally called ancilla, that can be in a 
superposition of the states |0) and 1 1). As shown in Fig. [3] the 
state of the ancilla controls the operation of the last Hadamard 
gate on the top line. In our implementation, we have chosen 
to include a preparation procedure for the state of the ancilla, 
as indicated in Fig. [3] 

For completeness and comparison with the event-by-event 
simulation data, we give the quantum-theoretical descrip- 
tion of this experiment in terms of the state \vu) = |v) (S> \u) 
where h,v = 0,1 label the basis states and \u) and |v) de- 
note the state of the ancilla and photon, respectively. The 
amplitudes at the input a = (aoO)«Oi)«io,aii) r and output 
b = (^oo,^Oi7^io,^ii) r of the experiment depicted in Fig. [3] 
are related by 



More explicitly we have [ 34] 
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(1) 



where a = 1 



Reading from right to left, the matrices in Eq. ([T) represent 
the action of a Hadamard operation on the ancilla, a phase 
shift (by a) operation on the ancilla, another Hadamard oper- 
ation on the ancilla, a Hadamard operation on the photon, a 
phase shift (by <p) operation on the photon, and a controlled 
(by the ancilla) Hadamard operation on the photon. Note that 
all these operations only affect the state, that is the wave func- 
tion, which describes the statistical properties of the whole 
system and cannot be interpreted as having causal effects on a 
particular particle without running into conceptual and logical 
problems |4|]. 

For the case at hand, «oo = 1 and all other a's are zero. Then 
it follows from Eq. ([TJ that the probability to detect a pair 
(photon, ancilla) in the state \vu) is given by p(v,u) = \b VM \ 2 . 
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Note that Eq. (f2|i is identical to the corresponding result for 
the standard delayed-choice experiment. 



SIMULATION MODEL 



The model presented in this paper builds on earlier work 19- 



el presentei 

lli|2j|23j|27iIol|3i in which we have demonstrated that it 
may be possible to simulate quantum phenomena on the level 
of individual events without invoking concepts of quantum 
theory. 

In our simulation approach, a messenger (representing the 
photon or the ancilla), carries a message (representing the 
phase) and is routed through the network and the various units 
that process the messages. 

We now explicitly describe our simulation model that is, 
we specify the message carried by the messengers, the algo- 
rithms that simulate the processing units and the data analysis 
procedure. 

Messenger. Particles carry a message represented by a 
two-dimensional unit vector y^ „ = (cos y/ J t „,siny/ / (-„) where 

„ refers to the phase of the photon. The subscript n > 
numbers the consecutive messages and k = 0, 1 labels the 
port of the beam splitter at which the message arrives. Ev- 
ery time, a messenger is created, the message is initialized to 

Hadamard gate. The key element of the event-by-event 
approach is a processing unit that is adaptive, that is it can 
learn from the messengers that arrive at its input ports lloT Hlll . 
The processing unit consists of an input stage called deter- 
ministic learning machine (DLM) 19j llOfl, a transformation 
stage, and an output stage. In experiments with single parti- 
cles, the input stage receives a message on either input port 
k = or k = 1 , but never on both ports simultaneously. The 
arrival of a message on port 0(1) corresponds to an event of 
type 0(1). The input events are represented by the vectors 
e„ = (1,0) or e„ = (0, 1) if the «th event occurred on port 
or 1, respectively. The DLM has two sets of internal regis- 
ters (Q; ,,,^^) and one internal vector x„ = (xo, H ,xi. n ), where 
x o.n+ x i,n = 1 an d x i,n > 0. These three two-dimensional 
vectors are labeled by the message number n because their 
content is updated every time the DLM receives a message. 
Thus, the DLM can only store 6 numbers, not more. Before 
the simulation starts we set xo = (*o,o,xi,o) = (<^>1 — 
where 2$, is a uniform pseudo-random number. In a simi- 
lar way, we use pseudo-random numbers to set (Cto,St o) 
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for k — 0, 1. Upon receiving the (n + l)th input event, the 
DLM performs the following steps: (1) it stores the mes- 
sage ytn+i = (cos % >n+ b sin in its internal register 
(Q-,n+i!^<r,n+i) an d (2) it updates its internal vector accord- 
ing to the rule 



Xi,n+\ = JXi.n + (i - 



(3) 



where < 7 < 1 is a parameter that controls the learning pro- 
cess. By construction .^o.n+i +%in+l = 1 and;t,„ + i > 0. 

The parameter 7 affects the time that the machine needs to 
adapt to a new situation, that is when the ratio of particles on 
paths and 1 changes. By reducing 7, the time to adapt de- 
creases but the accuracy with which the machine reproduces 
the ratio also decreases. In the limit that 7=0, the machine 
learns nothing: it simply echoes the last message that it re- 
ceived JsjlQii. If 7^ 1~, the machine learns slowly and ac- 
curately reproduces the ratio of particles that enter via path 
and 1. It is in this case that the machine can be used to 
reproduce, event-by-event, the interference patterns that are 
characteristic of quantum phenomena 

The transformation stage implements the specific function- 
ality of the unit, the Hadamard operation for the case at hand. 
It takes as input the data stored in the two internal registers 
(Ck, n +i,Sk,n+i) (k — 0,1) and in the internal vector x„ + i = 
C*0,n+i,-*i,n+i) an d constructs the four-dimensional vector 



^ Cl,H+l\Al,n+l +Co,„+ly/XQ. n+ l \ 

S\,n+1 y/ x l,n+l +So.n+ \^f x 0.n+\ 
\ S\.n+l^/Xl,n+\ — S(ln+1 A /*0,M+1 / 



(4) 



Rewriting this vector as a two-dimensional vector with 
complex-valued entries, it is easy to show that V corresponds 
to the matrix-vector multiplication in the quantum theoretical 
description of the Hadamard gate B23I1 . 

The vector V is then passed to the output stage which de- 
termines the output port through which the messenger leaves 
the unit. The output stage sends the message 



W0, I+ l = (V + V 1 )/(V2 + V2) 1 /2 7 

through output port if Wq ji+1 < ^ where < 



(5) 



< 1 is a 



uniform pseudo-random number. Otherwise, the output stage 
sends the message 



Wi„+i = (V 2 - 



•v 3 )/(v|- 



. V 2 ) l/2 j 



(6) 



through output port 1 . 

Controlled Hadamard gate. The event-based processor of 
this device is identical to the one of the Hadamard gate itself 
except that the vector V is computed according to Eq. dU if 
the control bit (called x) is 1 only. If the control bit is 0, V is 
given by 



/ Co,«+i\/ x o.«+i \ 

S0,n+ly/X(ln+l 
Cl,«+1 V^l.n+l 



(7) 



Phase gate. The unit that performs the phase shift by an 
angle (j) changes the message „ according to the rule 



yo.« 



yi,« <- 



1 
1 



COS0 

-sin0 



sin0 

COS0 



yi. n - 



(8) 



As a result the message is rotated by (j) if the particle traveled 
via path 1 . 

Simulation procedure. For each pair (ct,(p), N = 10000 
pairs of messengers (one for the photon, one for the ancilla) 
are sent through the network (see Fig. [3]l of processing units. 
A messenger that appears on an output line of the network, 
either exits via port or via port 1, never via both ports simul- 
taneously. With each pair of messengers that emerges from 
the network, the corresponding counter is incremented, that is 
no events are being discarded. In other words, we assume that 
the efficiency of the detectors is 100%. After all pairs have 
been processed, dividing the value of one of the counters by 
N yields the normalized frequency for observing a pair (pho- 
ton,ancilla) in the corresponding output ports. 



SIMULATION RESULTS 

In Fig. |U we show results of the event-based simulation 
of the quantum-controlled delayed-choice experiment for a 
fixed value (a = tt/3) of the parameter that determines the 
probability (sin 2 a) that the ancilla is in the state |1). As the 
solid lines in Fig. [4] are the predictions of quantum theory, see 
Eq. (0, it is clear that the event-based simulation reproduces 
the results of quantum theory for this particular value of a. 

In Fig. [5] we plot the difference between the event-based 
simulation results and the prediction of quantum theory, given 
by Eq. ©. The differences are on the 1% level, as it should be 
on the basis of standard statistical arguments. Therefore, we 
may conclude that the event-by-event approach reproduces the 
statistical distributions of quantum theory for the quantum- 
controlled delayed choice experiment. 



DISCUSSION 

Instead of discussing our event-by-event simulation ap- 
proach for optical phenomena in full generality, in this pa- 
per we have opted to explain in detail how the approach is 
applied to a specific example, a quantum-controlled delayed- 
choice experiment 113411 . We hope that this helps to understand 
the key feature of our approach, namely that it builds, one- 
by-one, the statistical distributions of quantum theory without 
knowing about the latter. 

The successful simulation of the quantum-controlled 
delayed-choice experiment 13411 adds one to the many exam- 
ples for which the event-by-event simulation method yields 
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<p (degrees) 

FIG. 4: The normalized frequency of observing a photon in path 
(squares) or 1 (circles) conditioned on the observation of the ancilla 
in path (open symbols) or 1 (closed symbols), for the case in which 
a = 7r/3. The solid lines are the prediction of quantum theory, see 
Eq. (|2j . The number of emitted and detected events per (f> is 10000. 
The DLM control parameter y = 0.99. 




FIG. 5: Difference A(a,<p) between the quantum theoretical result 
Eq. (O and the data obtained from an event-by-event simulation of 
the quantum-circuit shown in Fig. [3] The number of emitted and de- 
tected events per pair (a, <J>) is 10000. The DLM control parameter 
7= 0.99. The differences fluctuate on the 1% level. Open squares: 
Photon detected in path 0, ancilla detected in path 0; Closed squares: 
Photon detected in path 1, ancilla detected in path 0; Open circles: 
Photon detected in path 0, ancilla detected in path 1; Closed circles: 
Photon detected in path 1, ancilla detected in path 1. Lines connect- 
ing markers are guide to the eye only. 



the correct statistical distributions. Of course, the event- 
based approach, being free from concepts such as particle- 
wave duality, does not suffer from the conflicts with every- 
day logic that arise in the quantum-theoretical description of 
the delayed-choice experiment. In particular, there is no need 
to invoke the thought that in this experiment, the character of 
the photon need to be changed in the past. 

Finally, it should be noted that although the discrete-event 



algorithm can be given an interpretation as a realistic cause- 
and-effect description that is free of logical difficulties and 
reproduces the statistical results of quantum theory, at present 
the lack of relevant data make it impossible to decide whether 
or not such algorithms are realized by Nature. Only new, ded- 
icated experiments that provide information beyond the statis- 
tics can teach us more about this intriguing question. 
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